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PREFACE 

The main purpose of Tn No* 64T-1 is to serve as a guide to students 

using the Fortran cards for computer program "CMPRSOR" to reduce the 
data taken during a classroom test of the compressor test rig. This 
compressor test rig is described in USNPGS Technical Report No, 12, 

To reduce data taken during a compressor test rig run y it is 
necessary only to follow directions on pages 7, 8 and 9 of this technical 
note. Appendix A will also help toward an understanding of the program 
and its use. 

Appendix B gives the information necessary for slide rule calculations 

of efficiency data. 

In content, TN No 0 64T-1 is a compilation of information revelant 
to the compressor test rig not contained in USNPGS Technical Report No. 12, 
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1. SYMBOLS : 



A 

S 

T 

T 

V 

W 

W 

W 

W 



F 

c 



1 

2 

L 



Amps . 

Slip - RPM 
Temperature - op 

Temperature - C 

Volts 

Watts 

Watts (Left Hand Watt meter) 

Watts (Right Hand Watt meter) 
Watts Loss 
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2. INTRODUCTION : 

Technical Report No. 12 published in February 1955 by Professors 
Vavra and Gawain covers the construction and use of the compressor 
test rig and gives all information required for data reduction leading 
to values for flow and pressure coefficients cp and \|r. The derivation 
fo formulas for flow and pressure coefficients is also included in TR No. 12. 
The information given for slide-rule computation of flow and pressure co- 
efficients cp and ; is still pertinent; however, the computer program des- 
cribed and discussed in the above report applied to a computer no longer 
at the USNPGS. 

The following supplements Report No 12, in two respects: First* 

it shows the development of formulas for calculation of compressor 
efficiency by the use of a computer program and* second, it gives detailed 
information on the use of a Fortran program^ This program is for the 
reduction of all data taken during a compressor test rig run. 

3. COMPUTER PROGRAMMING OF MOTOR OUTPUT DATA : 

Before 1956 the compressor test rig was operated without power input 
instrumentation. In 1956 the Westinghouse Company calibrated the compressor 
drive motor and supplied loss curves. 

Net input to the compressor is equal to the input to the motor minus 
losses, since the compressor rotor is connected directly to the motor and 
the motor bearings carry all bearing loads. Input to the motor is the 
algebraic (since one wattmeter may have a negative reading) sum of readings 
from two wattmeters, the”two-wattmeter system’* for 3-phase current. Total 
losses may be compiled by using the loss curves. Compressor efficiency, 
then, is the theoretical power required by the compressor divided by the 
net input to the compressor. A sample calculation (for slide-rule computa- 
tion) of the motor output and curves from data supplied by the Westinghouse 
Company are given in Appendix B 

In order to use the curves as part of a computer program it was necessary 
to reduce them to formulas Two of the curves are straight lines. Another 
curve has only three points. It was assumed that this curve is parabolic 
so the formula for the curve was set up as: 

y = A + BX + CX 2 
-A- 
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or, specifically: 

2 

Watts Loss = A + BX Volts + CX (Volts) 

Substituting values at the three points given, three equations can be 
written from which the unknowns A, B, and C can be computed^ 

The remaining curve has seven data points. With the exception of 
one point, these points are on a parabolic curve determined by any 
three of the six points. The points were checked and found to lie 
on this theoretical curve to an accuracy greater than the accuracy 
to which the curve drawn through the points can be read. 

For values above 12 amperes the constants of another parabolic 
curve were computed, using the 12 ampere point, the 14 ampere point and 
a point obtained by graphical interpolation between the 12 ampere point 
and the 14 ampere point. 

The computer program calls for the computer to choose which formula 
to use, depending on whether the input is above or below 12 amperes. With 
Fortran a choice is called for by the "IF" statement, 

A program for motor inpuf and loss, and thus for compressor efficiency, 
was first written to supplement the original compressor test rig program 
written for the CDC 102. Reduction of efficiency data is now part of 
a complete data reduction program for the compressor written in Fortran. 

For slide rule computations: 

Efficiency * Theoretical Power of Compressor /Input to Compressor 

Flow x Pressure x A Constant/Electrical Input - (Motor Losses + 

t r> 

"T " is the power loss due to the resisting torque of the water seal 
R 

mechanism . 

Input to the motor is measured by two wattmeters 
W = W x + W^ 

At low power input the right hand wattmeter (W.>) may be nega tive ; 
if this is indicated, two leads are reversed and the reading is recorded 
as negative. 

For slide rule computations the meter readings must be multiplied 
by a factor to give true volt, ampere, and watt values. For use in 
compiling data fo<r the Fortran program "CMPRSOR" meter readings are 
recorded exactly as read. 
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The ampere meter has multiple scales. The correct amperage is 
obtained by reading the second scale from the bottom and multiplying by 
10. For slide rule computations an average istaken of the readings 
from the three meters. 



Multiplying Factors 

Amperes (Second Scale from Bottom) Multiply by 10 
Volts " " 20 

Watts " " 40 

For compiling motor losses reference may be made to the loss curves 
and sample calculations in Appendix A. 

As an alternate method of computing motor losses for slide rule 
calculation of compressor efficiency the loss curve formulas may be 
used . 

The combined friction, windage stray and core loss is a function 
of the voltage." 

WI = -21(k0 + 0.8929 V - 0.0001047 V 2 

L 

Stator copper losses are a function of the amperage. 

W = 2.73 A 2 when A < 12 
b 

W L = 4484 - 783A + 33. OA 2 
when A > 12 

This loss is modified by a factor which is a function of motor winding 
temperature . 

Factor - 0.904 + 0.00385 T c 

Rotor copper losses are an indirect function of the slip. To get 
these losses the gross watts input minus the stator losses is multiplied 
by a factor which is a function of slip. 

Factor = 0.000556 x Slip 

When the compressor is set up for taking rotor blade pressures there 
is an added loss, which is not a motor loss. This is the loss due to the 
resisting torque of the water seal mechanism. When rotor blade pressures 
are not taken this loss goes on the data card as zero. 

For computing purposes the torque loss can be converted to equivalent 
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= Weight on Scales x Lever Arm x RPS x 2 TT x A Consta^ 
Lever Arm = 6" = 0.5' 

RPS = RPM/ 60 
2tt - 6.2832 

Constant = 1.356 (watts/f t-lbs/sec) 

W = Wt. x RPM x 0.07095 

L 

4. DATA REDUCTION USING FORTRAN : 

Program CMPRSOR was written to reduce data from the compressor 
test rig. A set of program cards is available. 

With the facilities provided by the school it is necessary only to 
know how to fill out Fortran data forms so the data can be card~punched 

There are 21 items to be entered as data. The data cards must be 
in the following order: 

First: The one card with 7 items of information. 

Second: 1 to 4 cards with "RC(J)" values 

Third: The two cards with the 12 data item- required for flow and 

efficiency computations at the first flow rate. 

Fourth: 1 to 4 cards with manometer board readings at the first flow 

rate . 

Repeat "Third" and "Fourth" for "N" number of flow rates. 

One line on the Fortran form equals one card, except that the Fortran 
form has only 12 spaces and the card has 80. Where all 80 spaces on a 
card are used, two lines on the Eprtran form are required For this prq'f 
gram 70 spaces (f<pr 7 items) can be on one line and the last data item 
can be entered on the first 8 spaces on the second line. Eighty card 
spaces are used 'for entering manometer readings and the n RC(J)" values. 

All data entries are entered in 10 spaces per data entry item -- thur 





1 


2 


3 


4 


5 


6 


7 8 9 10 


BA = 30.06 is 


3 


0 


. 


0 


6 






RPM = 1795 is 


1 


7 


9 


5 


• 






PZ = 2.20846 is 


2 


• 


2 


0 


8 


4 


6 



Floating point numbers may be placed anywhere in the assigned space. 
They are shown here in the first spaces. T$is system of palceraent: 
is the most convenient for the card punch operators. A decimal point 



- 7 - 



TN No. 64T-1 



must always be included for floating point numbers. 

The data is coded for ease of memory 'thus : : 

BA = barometer reading 
EP = e (epsilon) 

Six items to be entered on the first card: 

M = number of pressure readings for each flow rate. 

N = number of flow rates. 

NOTE: M and N are integers and must therefore be entered at the 

right-hand end of the space rpovided. IF M = 10 for 10 
pressure readings, enter it as: 

123456789 10 

1 0 

There is no decimal $oint. 

BA = barometric pressure., in inches of mecury. 

PZ - <P 0 — See p. 74 of ral £). Valu e. depend on orifice 
EP = e --- See p. 74 of TR12 diameter 

TARE 1 = Tare of micromanometer for P^ and Ap . 

TARE 2 = Tare of micromanometer for A 42. 

Next cards 1 to 4: 

RC(J) Correction for centrifugal force of air in tubes for rotor pressures 
Coded "Radius Correction". Find values on p. 75 of TR13U Enter "M" values 
For pressures not on rotor blades enter RC(J) as zero. One blank card 
enters 8 zero values. (Blank cards used here should not be confused with 
the one blank card between program cards and data cards.) 

For each flow value the following items are entered on two successive 
cards, the first 7 on the first card, the second 6 on the second card, 

1. W1 = Watts reading on left-hand meter 

Note: The meter reading (A, V, W) are entered as read. 

The program multiplies by the correct factor. 

2. W2 = Watts on right-hand meter. This may be negative. 

3. A = Amperes. Enter average of 3 readings . 

4. V = Voltmeter reading. 
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5. TW = Temperature of winding f - °C . 

6. WT = Weight in pounds recorded on scales. This is for measurements of 

rotor pressures only. If other pressures enter WT as zero. 

7. RPM = RPM 

1. PI S P 1 t 

2. T1 = Temperature in discharge pipe -- °F. 

3. PA = ha = Atmospheric reading on manometer board in centimeters of 

water . 

4. DP = £p’ = Inches of water. 

5. A42 = Plenum pressure reading at point A42 . Read from micromanometer 

±n inches of water. 

6. TA = Temperature of air entering compressor -- °F» 

"M" manometer readings are entered for n N n flow rates. PP(J) on computer 
program. Centimeters of water. 

For each flow rpte: 

1-8 readings - 1 card 

9 - 16 readings - 2 cards 

17 - 24 readings - 3 cards 

25 - 32 readings - 4 cards 

The table of results is self-explanatory. The pressure (\|f) values 
are printed according to the order .in which corresponding manometer reading 
pp(j>, were entered as data. 
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PROGRAM CMPRSOR 



1000 

2000 

3000 



20 



30 



40 



50 

10 

100 

300 

500 

600 

700 

800 

9000 



PROGRAM CMPRSOR 

DIMENSION PSI (96), KC(96), PP(96) 

READ’ 100, M, N, BA, PZ, EP, TAREl, TARE 2 
READ 500, (RC( J) , J “ 1, M) 

EBA - BA* (13.55) 

DO 10 I - 1, N 

READ 500, Wl, W2, A, ¥, TW, WT, RPM 

READ 500, PI, Tl, PA, DP, A42, TA 

TTA “ TA+459 . 4 

PP1 * P1“TARE 1 + EBA 

TT1 s Tl + 459.4 

DP * DP-TARE 

A2 ® FZ*(1.-EP*DP/PP1) 

A3 = 1800. /RPM * PP1/BBA % TTA/TTl 
A4 “ SQRTF(TT1/520.*DP/PP1) 

PHI - A2 * AS * A4 

OUT *» (A42 •fe- TARE 2) * PHI * RPM 

V = V * 20. 

A 53 A * 10. 

WA “ 40. * (Wl + W2) 

A2 = -210. +.8929 * V- . 0001047. * V * V 
IF 33 (A- 12.) 1000, 1000, 2000 
A3 -2.73* A *A 
GO s TO 3000 

A± s 4484.-737. * A+33. * A * A 
A3 - A3 * (.904 + .00385 * TW) 

A4 “ .000556 * (1800. -RPM) * (WA“A3) 

A5 - WT * RPM * .07095 

ET£.« (OUT/ (WA- (A2 + A3 + A4 + A5)>) * 151.196 
READ 500, (PP(J) 

DO 20 J = 1, M 

PSI (J) “ ( (PA"PP(J) ) * TTA (BA*RPM*RPM) * 64. 10.2) + RC.(J) 

WRITE OUTPUT TAPE 2,700, PHI 

WRITE OUTPUT TAPE 2,800, ETA 

WRITE OUTPOT TAPE 2,300 

WRITE OUTPUT TAPE 2,900 

WRITE OUTPUT TAPE 2,600, (PS1(J), J= 1,8) 

IF (M-8) 10, 10, 30 
WRITE OUTPUT TAPE 2,1100 

WRITE OUTPUT TAPE 2, 600, (PSI(J), J=9 , 16) 

IF (M-16) 10, 10, 40 

WRITE OUTPUT TAPE 2,1200 

WRITE OUTPUT TAPE2 ,600 , (PSI(J) , J= 17,24) 

IF (M-24) 10, 10, 50 
WRITE OUTPUT TAPE 2, 1300 

WRITE OUTPUT TAPE 2, 600, (PSI(J), J= 25,32) 

CONTINUE 
El® FILE 2 

FORMAT (2110, 5F10.5) 

FORMAT (27H PSI (PRESSURE COEFFICIENT)//) 

FORMAT (8F 10.5) 

FORMAT (8F9.5) 

FORMAT (/// /26H PHI (FLOW COEFFICIENT) = F6.5 //) 

FORMAT (2 OH ETA (EFFICIENCY) = F7.3.9H PER CENT//) 

FORMAT (68H 1 2 3 4 5 

7 8 ) 



6 
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PROGRAM CMPRSOR (Cont r d) 



! 



11000 

1 


FORMAT 

15 


(69H 


9 

16) 


IP 


11 


12 


13 


14 


12000 

1 


FORMAT 

23 


(69K 


17 

24) 


18 


19 


20 


21 


22 


13000 

1 


FORMAT 

31 


(69H 


25 

32) 


26 


27 


28 


29 


20 



END 

END 
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SAMPLE CALCULATION 
Horsepower Output of Motor 



Measured Test Data; 



(1) Input^oltage 

(2) Input Amperes 

(3) Input Watts 

(4) Winding Temp. 

(5) RPM 



(ave. of 3 legs) 
(ave of 3 legs) 

(°C) 



Motor Losses 



from Curves; 




- 2 300 

= 10 

= 31760 
= 20 

= 1795 



( 6 ) 

( 7 ) 

( 8 ) 

(9) 

( 10 ) 



( 11 ) 

( 12 ) 

(13) 



Combined Friction - Windage etc. 

Stator Copper Loss at 25°C = 273 

Temperature Correction Factor * .981 

Stator Copper Loss at Winding Temp. 273 x .981 

Slip Correction Factor .00834 

Input Watts 31760 

Minus Stator Copper Loss 268 

Input Factor 31492 

Rotor Copper Loss 

Input Factor 3.492 (11) x Slip Corr. .00834 (10) 
Total Motor Losses (6) + (9) + (12) 



Motor Output; 

(14) Input 31760 (3) - Losses 1820 (13) 

(15) Watts Output (14) x .001341 = HP 



- 1290 



= 268 



» 262 
1820 



= 29940 
- 40.2 
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Notes on the U$e of the Compressor Test Rig for a laboratory test 

1* The oil in the bearings should be checked when starting for the first 
time after, say, a week of no running. For shorter periods of time, loss 
of oil would be indicated by pools of oil underneath the motor. The bearings 
run cool so any bearing heat would indicate trouble. The motor should 
take from three to five minutes to come to a complete stop. 

2 There has been no bearing trouble of any kind. However, there 
is an indication that bearing friction is a considerable motor loss. When 
the electrical instrumentation was first installed it was found that the 
input to the motor decreased during the first six hours of continuous running. 
One curve shows 4 decrease of 8.37, in the first four hours. This was 
analysed as due to lowered bearing oil viscosity because of increased oil 
temperature. If possible the compressor should be run for two hours just 
before a laboratory clasps :i^ held. 

2. The motor must be started with the switches to the low-resistance 
shunts, that are across the instrument leads , closed. If they are left 

open the high starting amperage can ruin an instrument . There is also a switch 
on each of the meters that must be closed in order to read the meter. 

These switches may be locked open or they may be held closed for the time that 
the meter is being read. 

3. ^he sharp-edge orifice in the discharge pipe used for flow measurements 
has an inside diameter of 13.335 n o If the orifice is* as before, taken out 
and the edge sharpened by a machining operation , the inside 'ftj.ameter must be 
measured accurately and the value recorded for flow measurement calculations. 

4. There are several possibilities for pressure leaks in the lines from the 

orifice flange to the micromanometer: first, the connections at the orifice 

flanges may leak and, second, there is a possibility of pressure leaks at 

the quick-disconnect arrangement (the brass block and rubber cork rig) and, 
finally, the connections at, or even inside, the micromanometer, afap leak. 

For a rough check on the possibility of leaks it is noted that, at full flow, 

Ap is about 14" of water for 30° conterflow setting of the inlet vanes and 
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about 10" for zero setting (vanes radial) • T his is for the usual 
40° setting of the vanes at the rotor discharge and for the 13.535 n flow 
orifice, which is the one used. The 13.535 n orifice was itiginally 
13.500" I.D. The 11.5" and 15.5" orifices are not used. 

5. For a test the following items must be at hand 

A thermometer 

Two micromanometers 

Four clip boards 

Blank data sheets „ There are kept in Room 123, Bldg. 234. 

RPM counter 

USNPGS Technical Report No. 12 

Several copies of this publication 

Ozalid print of cross-section of test rig, drawing No* 665o 

Ozalid print of cross-section of water seal mechanism, drawing No. 542B* 

One micromanometer is used for flow measurement. 

One micromanometer is connected to pressure tap A42 for the plenum 
pressure. 

t Pressure lines will be connected to a 40- tube manometer board as 
outlined by the professor conducting the class. 

The ozalid prints are to aid the professor holding the class in explaining 
the rig and its use. 

6. When taking rotor vane pressures great care must be taken in entering 
the water used for isloating the pressure chatobers. If the water is let in 
at too high a head some of it may get into one or more ot the small- 
diameter pressure lines leading to the tap locations. When this happens 

a reading at the tap affected cannot be taken until the rotor is removed 
and the line cleared. The can that is raised fo force water into the water 
seal mechanism is raised Slotf&y two or three feet above the center- line 
of the compressor and lowered as soon as all chambers are isolated. A 
manometer board pressure reading not the same as that from either of the 
two adjacent isolation chamber readings, and constat, is an indication 
that sufficient water has entered to isolate the chamber connected to that 
particular tube. More than sufficient water may cause water to be trapped in 
the pressure lines from the rotor tapeldcations . 
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